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We reconstruct dark energy properties from two complementary supernova datasets - the newly 
released Gold+HST sample and SNLS. The results obtained are consistent with standard ACDM 
model within 2a error bars although the Gold+HST data favour evolving dark energy slightly more 
than SNLS. Using complementary data from baryon acoustic oscillations and the cosmic microwave 
background to constrain dark energy, we find that our results in this case are strongly dependent on 
the present value of the matter density fiom. Consequently, no firm conclusions regarding constancy 
or variability of dark energy density can be drawn from these data alone unless the value of Qo m 
is known to an accuracy of a few percent. However, possible variability is significantly restricted if 
this data is used in conjunction with supernova data. 



INTRODUCTION 



There is a growing consensus in cosmology that the Universe is currently accelerating. Perhaps the simplest 
explanation of this property is the presence of a positive cosmological constant A. Although A appears to explain 
all current observations satisfactorily, to do so its value must necessarily be very small A/8ttG ~ 10~ 47 GeV 4 . So, 
it represents a new small constant of nature in addition to those known from elementary particle physics. However, 
t-H | since it is not known at present how to derive it from these other small constants and it is even unclear if it should 
be exactly constant, other phenomenological explanations for cosmic acceleration have been suggested. Collectively 
called dark energy (DE) models, they are based either on the introduction of new physical fields (quintessence and 
phantom models, the Chaplygin gas, etc.), or on geometrical approaches which attempt to generate acceleration by 
means of a change in the laws of gravity and, therefore, the geometry of space-time Scalar-tensor gravity, R+f(R) 
■ gravity and higher dimensional 'Braneworld' models are prominent members of this second category. 

The growing number of DE models has inspired a complementary approach whose aim is to reconstruct properties 
O ,■ of DE directly from observations in a quasi-model independent manner, see the recent review Q and the extensive 
list of references therein. The aim of this paper is to investigate what new insights can be obtained about DE using 
the most recent data, and to see whether or not these data strengthen previously obtained results on the closeness of 
DE to A. We study two different supernova samples - the newly released Gold+HST sample [1| and SNLS data [1] - 
to see what constraints follow from them on possible evolution of DE. We also investigate the possibility of extracting 
information about the nature of DE from datasets other than the supernova data. To this end, we consider what 
follows from the value of the R parameter characterizing acoustic peaks in the angular power spectrum of the cosmic 
microwave background (CMB) @, 0| and from the SDSS result on the baryon acoustic oscillation peak (BAO) 



DATA AND METHODOLOGY 



In this paper, we shall compare the reconstruction results for different recent sets of observations. We briefly 
summarize each of the data sets which we shall use before proceeding to give the results of our comparison. 

Type la supernova data is the strongest evidence for DE in current cosmology. The data is in the form of [M),i(zi) = 
wis — M = 5logd,L(zi) + 25 with 

d L (z) = (i + z) [ 4h> W 

Jo H ( z ) 

where H(z) = Hoh(z) is defined in ([8]) for the ansatz used. 

We use the following two SNe datasets available to us at present. 

1. The Gold+HST data set : As recently as 2003, the entire supernova dataset from the two different surveys 
- Supernova Cosmology Project (SCP) and High z Supernova Search Team (HZT), along with low redshift 
supernovae from Calan-Tololo Supernova Search (CTSS) comprised of a meager 92 supernovae d, 0, [l(| , with 
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very few at high redshifts, z > 0.7. The method of data reduction for the different teams was also somewhat 
different, so that it was not possible to use the supernovae from the two datasets concurrently. The picture 
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changed somewhat dramatically during 2003-2004, when a set of papers from both these teams 
presented a joint dataset of 194 supernovae which used the same data reduction method. This new data 
resulted in doubling the dataset at z > 0.7. Not all these supernovae could be identified beyond doubt as 
Type la supernovae however, in many cases complete spectral data was not available. In early 2004, Riess et 
al. 1J] reanalyzed the data with somewhat more rigorous standards, excluding several supernovae for uncertain 
classification or inaccurate colour measurements. They also added 14 new high redshift supernovae observed 
by the Hubble Space Telescope (HST) to this sample. This resulted in a sample known as the 'Gold' dataset. 
The latest publication from the same team Q adds a further 10 SNe from HST to the dataset, and excludes 
data below cz = 7000fcm/ s (z — 0.0233) to avoid the influence of a possible local "Hubble Bubble". The final 
sample now comprises of 135 supernovae (the furthest being at redshift z = 1.755) - we call this sample the 
'Gold+HST' dataset and use it as our first SNe sample. 



2. The Supernova Legacy Survey SNe data set (SNLS) : The SuperNova Legacy Survey [15j is an ongoing 5-year 
project which is expected to yield more than 700 spectroscopically confirmed supernovae below redshift of one. 
The first two-year results from this survey [3] have provided us with 71 new supernovae below z = 1. We shall 
use these 71 SNe together with the already available low- 2 supernova data, i.e. a total of 115 SNe, as our second 
supernova sample. 

In the SNLS and Gold datasets have also been merged together for some of the results. Since the standardization 
techniques of the two teams are somewhat different, as of now, we consider the two samples separately till a better 
handle on the systematics is obtained. 

Although supernova data forms the major observational proof of dark energy, we may try to obtain independent 
information on the nature of dark energy from the other observational results. Here we look at the following two 
datasets. 

1. Observations of the cosmic microwave background provide us with very accurate measurements, which may be 
used to gain insight about dark energy [f|. We may use the WMAP 3-year results to get 0] 
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= 1.70 ±0.03 
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where h(z) = H(z)/Ho is defined in (JSJ) for the ansatz used. For calculating this quantity, we use H,t,h 2 = 0.024 
and fiom's 2 = 0.14 ± 0.02. To calculate zi s we use a fitting function given in (l6j : 



z ls = 1048[1 + 0.00124(lVi 2 )- - 738 ][l + gi (n Qm h 2 ) 92 ] 
where the quantities g\ , gi are defined as 

9l = 0.078(f^ 2 r - 238 [l + 39.5(fM 2 ) - 763 r\ 



.92 = 0.56[l + 21.1(fW) 
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2. The Baryon Acoustic Oscillation Peak : A remarkable confirmation of the standard big bang cosmology has 
been the recent detection of a peak in the correlation function of luminous red galaxies in the Sloan Digital Sky 
Survey This peak, which is predicted to arise precisely at the measured scale of 100 h _1 Mpc due to acoustic 
oscillations in the primordial baryon-photon plasma prior to recombination, can provide a 'standard ruler' with 
which to test dark energy models. Specifically, we shall use the value [7j 



A = 



Om 



h( Zl y/3 



i 



Zl dz_ 



-i 2/3 



0.469 



(—) 
V0.98/ 



-0.35 



±0.017 



(6) 



where h(z) — H(z)/Hq is defined in © for the ansatz used, and z\ = 0.35 is the redshift at which the 
acoustic scale has been measured. The 3-year WMAP results yield n — 0.95 for the spectral index of density 
perturbations. 



We note here that the quantity A measured from the baryon acoustic oscillations appears to be independent of 
the value of h. However, the quantity R obtained from CMB measurements depends on flo m h 2 , and therefore on the 
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Hubble parameter. Therefore, when using this quantity, we require to either assign some value to h or to marginalize 
over £lo m h 2 . 

We use x 2 minimization on the different datasets with 

2 itt \ XT [yfit.»(^;#o,^Om,Pj) - Vi] 2 ,„s 

x (.■ttoiPj) = 2^, 72 ■ (') 

i * 

Here, yi is the data at redshift of Zi and Oi is the uncertainty in the individual y.;, and pj are the parameters f2om, Ai, A 2 



for the ansatz [17 1 

#, 2 



if 2 (z) 

fe 2 (z) = 2 = n 0m a: 3 + A + Aix + A 2 x 2 , x = 1 + z (8) 



with A = 1 — Slom — Ai — A 2 for a flat universe. Note that this ansatz for the Hubble parameter has been tested 
and found to give accurate results for a variety of cosmological models TH 1||. Its main features are: (i) it accounts 



for a matter dominated regime at high redshift, (ii) it is exact for DE being the mixture of a cosmological constant 
and non-relativistic domain walls and cosmic strings (in other words, for DE consisting of 3 components with with 
w = — 1, w = — 1/3 and w = —2/3), (iii) it can satisfactorily emulate DE with an evolving equation of state including 
Braneworld models and phantom models admitting regions with w < — 1 (moreover, the DE effective energy density 
Pde may even become negative in such models). 

RESULTS 

To test the two SNe datasets, we first reconstruct a ACDM model with curvature, using the formula 

d L (z) = ( f [(1 + + n m z>) - z'(2 + z')n A ] - 1/2 dz>) , (9) 

vM V Jo J 

where 

S! m + Cl\ > 1 => S(x) = sin(x), k = 1 — fi m — fl\ 
f2 m + < 1 =>• S(x) — smh(x), k = 1 — Q m — 51a 
O m + 51a = 1 => = x i k = 1 . 

In figure [1] we show results for a ACDM model for both the Gold+HST and SNLS datasets. We see that for the 
Gold+HST data, the dataset prefers a closed universe somewhat over the flat model. The SNLS best-fit, on the other 
hand, is very close to the flat model. 

We next study the two datasets in more detail using ansatz ([5]). In figure we show the results for the Gold+HST 
SNe dataset. We marginalize over flo m = 0.28 ± 0.03 which is the currently accepted value for matter density from 
SDSS and show the 2a confidence levels in Ai — A 2 in panel (a). The logarithmic variation of pue = A0 + A1X+A2X 2 
is shown in panel (b). The variation of the equation of state of DE 

. f/^';^- 3 ' (10 , 

1 - (Ho/Hyil 0m x A 

is shown in panel (c), and the deceleration parameter 

, 9 H'(x) , s 

q = -a aR 2 = — x - 1 , x = l + z 11 

is plotted in panel (d). We find that the best-fit is somewhat away from ACDM, and favours an evolving model of DE, 
although ACDM is still within 2a of the best-fit. From the deceleration parameter we find the value of the acceleration 
epoch, which is the redshift at which the universe started accelerating. For the Gold data this transition between 
deceleration and acceleration occurs between z a = 0.32 — 0.48 at 2a. The earlier Gold dataset of [141) favoured an 
evolving DE model starting from w < — 1 at present and evolving to w ~ at z ~ 1 (as shown in [20| and numerous 
subsequent papers), the addition of the 10 new HST data points do not appear to change that result. These results 
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also agree with the analysis shown in figure 12 of 3], as well as results from the non-parametric approach advocated 



Figure [3] shows the results for the same reconstruction for the SNLS dataset. Once again we marginalize over the 
matter density. The results in this case are much closer to ACDM, the DE equation of state shows very mild evolution. 
The acceleration epoch is in the range z a = 0.41 — 1.0. 

We note here that although we have marginalized over a particular matter density, results from the supernova data 
are only weakly dependent on Slom ■ Within the reasonable range f^om = 0.22 — 0.34 (the 2tr limits on the current 
observations for f^om), the results do not change significantly. 

The ansatz ([8]) is useful not only for reconstructing the DE density but also for determining an important related 
quantity - the w-probe [22] ]. As its name suggests, the w-probe provides us with important insights about the equation 
of state. However, unlike w(z) which is determined by differentiating H(z), the w-probe is related to the difference 
in the value of the Hubble parameter and is defined as follows 

1 , - 1 /* h , f m,, . n 1 A ln(p DE /poc) 

1 + W = Aln(l + z) J [1 + W[z)]d ln(1 + Z)= 3 Aln(l + z) ' (12) 

Here A denotes the total change in a variable between integration limits. It is important to note that the w-probe, 
when studied over different redshift bins, can give us evidence for the variation of the equation of state from the 



first derivative of the data alone, i.e. from the DE density. (As demonstrated in [22j, w is much less sensitive to 
uncertainties in the value of the matter density than w(z).) Table U shows the value of the w-probe in three redshift 
bins (approximately equal in log(l + z)) for both Gold+HST and SNLS datasets. We see that for the gold dataset 
w is close to —1 in the low redshift bin, but in the two higher redshift bins w is closer to zero. Thus the Gold+HST 
dataset seems to provide support for the evolution of DE. On the other hand, for the SNLS data, for both redshift 
bins (upto redshift of unity), w is close to —1, showing that this dataset favours the cosmological constant more than 
Gold+HST. 

In a way, the SNLS results contradict the results obtained for the Gold dataset. This fact had been noted for an 
earlier Gold dataset in [23| ■ However, the difference is not very large, at 2a the results from both datasets are consistent 
with each other. We should also keep in mind the fact that the light-curve standardization of the two datasets is done 
using different methods, MLCS2k2 for Gold+HST and SALT for SNLS. As shown in 0], different standardization 
techniques may lead to differences of upto 0.16 magnitude in the data. The discrepancy in cosmological results may 
therefore be attributed partly to the different standardization techniques. There are also possible effects from other 
sources in the data, such as systematic noise and K-correction. Therefore it seems that using just the two current 
SNe datasets, a firm conclusion cannot yet be reached about the nature of DE. We may conclude that while the 
cosmological constant is more or less consistent with both datasets, evolving DE is still not ruled out. 

Next we examine our results for the CMB+BAO data. In figure 2] we show results for this dataset for which the 
observations are at redshifts z = 0.35 and z = 1089. We present our results in panels (b), (c) and (d) over the redshift 
range of z = — 1.7 in order to be able to compare with our earlier results obtained using supernova data. The results 
for a matter density marginalization of f^om = 0.28 + 0.03 are reasonably consistent with ACDM, the confidence levels 
for A\ — A2 are somewhat tighter than those seen for the SNe results (Comparing panel (a) of figures [2j [3J S|) . However, 
unlike the SNe data, the results from CMB+BAO seem to depend much more strongly on the chosen value of the 
matter density. In figure [5] we show the difference in the results for three different marginalizations of the matter 
density. In the left hand panels, we marginalize over f2om = 0.25 ± 0.03, and the data then favours evolving DE 
which crosses the so-called phantom divide at wo = — 1. In the middle panels, we marginalize over £lo m = 0.28 + 0.03, 
the result is a mildly evolving DE model consistent with ACDM. The panels to the right show results marginalized 
over a higher value of the matter density Slo m = 0.31 ± 0.03, which suggests a DE model with practically constant 
W > — 1. These results are supported by an analysis of w-probe shown in table ITT1 which demonstrates that the 
redshift dependence of the mean equation of state w is sensitive to the value of Jlom- Thus, unlike the SNe data, even 
a very slight difference in fio m can lead to a significant difference in the degree of evolution of DE for the CMB+BAO 
dataset. Note that the generic trend of fl m dependence is similar to that recently found in [StJ where all data (the old 
Gold+HST sample, the SNLS dataset, WMAP, BAO and others) were analyzed using the Chevallier-Polarski-Linder 



paramctrization of the DE equation of state |28|, |29| . 

Finally we combine the SNe and CMB+BAO data. We use the marginalization f2o m = 0.28 ± 0.03 for the entire 
dataset. The results are shown for Gold+HST data in figure[6] We see that the parameter space is strongly constrained 
by the three sets of data, the cosmological constant is consistent with the data, so are mildly evolving DE models. 
The results from SNLS are shown in figure [7l We see that the results from both sets of SNe data are strikingly 
similar, this is because the results are very strongly constrained by the CMB+BAO data. Whether the CMB+BAO 
data should be used in conjunction with the SNe data, given the strong matter density dependence of the former, is 
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however rather questionable at this point. 

CONCLUSIONS 

In this paper, we have explored the reconstruction of DE using the most recent and complementary datasets 
available to us at present. First, we reconstructed DE using two supernova datasets- the recently released Gold+HST 
dataset and the SNLS two-year dataset. We find that the results for the two are slightly inconsistent. The Gold+HST 
dataset appears to favour an evolving model of DE with wq < —1 at present and w ~ at z ~ 1 over ACDM for 
its best-fit, however ACDM is still consistent at 2cr. For SNLS data, the results are more consistent with ACDM. It 
should be emphasized here that we are speaking about DE properties averaged over a redshift interval ~ 0.4 or more, 
see e.g. Tables I, II. Note also that, as follows from these tables and Figs. 2-7, significantly phantom behaviour of DE 
(w < —1-1), if it occurs at all, is possible for moderately low redshifts z < 0.4 only and is favoured by the Gold+HST 
dataset, but not by the SNLS data. 

This discrepancy may be a result of the different light-curve standardization techniques used by the two teams or 
due to systematic bias. Also, this could be due to the fact that the low redshift and high redshift SNe are obtained 
by different surveys. Future surveys such as SNAP [24(, JEDI [25| and DUNE [26| should have a better handle on 
the systematic errors and such discrepancies should disappear. 

We also attempt to obtain information on DE using datasets complementary to the supernova sample. To this end, 
we use information on the shift parameter R from the 3-year WMAP data and the quantity A for the baryon acoustic 
oscillation peak from SDSS. The results are consistent with ACDM but do not rule out weakly time dependent DE. 
These datasets appear to be rather sensitive to the value chosen for the present matter density f^om- As a result, it is 
difficult to reach firm conclusion on the nature of DE from these data until strong model-independent constraints on 
f2om are obtained. Thus, even with the most recent data, the fundamental question if DE reduces to A or not, still 
remains open. On the other hand, possible deviation of DE properties from those of A is gradually becoming more 
and more restricted. 

While this work was being finalised the paper [30| appeared containing results which are in broad agreement with 
ours in areas of overlap. 
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TABLE I: The reconstructed ii>-probe for the two different SNe datasets, Gold+HST and SNLS, using fi m = 0.28 ± 0.03 

W 

Dataset Az = - 0.414 Az = 0.414 - 1 Az = 1 - 1.755 



Gold+HST -1.160lH 8 7o -0.226t°;2iS 0.2681 ,; ,!? 



SNLS -1.0371°;°?! -0.985±g-.g| 



TABLE IL The reconstructed u>-probe for the CMB+BAO dataset using three different marginalizations over Qo n 

W 

n 0m Az = - 0.414 Az = 0.414 - 1 Az = 1 - 1.755 
0.25 ±0.03 -1.1951°;°?° -0.9481°;™ -0.493t°;°£ 
0.28 ±0.03 -0.9511°;°™ -0.779t°;|£ -0.6291°,;^ 
0.31 ±0.03 -0.772l°;°?9 -0.8191°,;^ -0.8521°,;^ 




FIG. 1: 2a confidence levels are shown for the Gold+HST dataset (left panel) and SNLS dataset (right panel) in the fiom,f2A 
plane, with the dotted line showing the flat universe. 
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FIG. 2: 2a confidence levels for the Gold dataset using fiom = 0.28 ± 0.03. The upper left hand panel shows the confidence 
levels in Ai — A2, with the black dot representing ACDM. The upper right hand panel shows the logarithmic 2a variation of 
the DE density in terms of redshift. The dashed line represents ACDM. The lower left and right hand panels represent the 
variation of the equation of state and deceleration parameter respectively. The dashed lines in both panels represent ACDM. 
The thick solid line in the lower right hand panel shows the acceleration epoch, i.e. the redshift at which the universe started 
accelerating. 
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FIG. 3: 2a confidence levels for the SNLS dataset using fiom = 0.28 ± 0.03. The upper left hand panel shows the confidence 
levels in Ai — A2, with the black dot representing ACDM. The upper right hand panel shows the logarithmic 2a variation of 
the DE density in terms of redshift. The dashed line represents ACDM. The lower left and right hand panels represent the 
variation of the equation of state and deceleration parameter respectively. The dashed lines in both panels represent ACDM. 
The thick solid line in the lower right hand panel shows the acceleration epoch, i.e. the redshift at which the universe started 
accelerating. Results are shown upto z = 1.01 
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FIG. 4: 2a confidence levels for the CMB+BAO dataset using f2o m = 0.28 ± 0.03. The upper left hand panel shows the 
confidence levels in Ai — A2, with the black dot representing ACDM. The upper right hand panel shows the logarithmic 2a 
variation of the DE density in terms of redshift. The dashed line represents ACDM. The lower left and right hand panels 
represent the variation of the equation of state and deceleration parameter respectively. The dashed lines in both panels 
represent ACDM. The thick solid line in the lower right hand panel shows the acceleration epoch, i.e. the redshift at which the 
universe started accelerating. 
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FIG. 5: 2a confidence levels for CMB+BAO using different marginalizations over fiom. The upper and lower left hand panels 
respectively show the confidence levels in A\ — A2 and the 2a variation of w over redshift for flom = 0.25 ± 0.03. The middle 
panels show the same for Qom = 0.28 ± 0.03 and the right-hand panels are results marginalized over fiom = 0.31 ± 0.03. The 
black dots in the upper panels and the dashed lines in the lower panels represent ACDM. 
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FIG. 6: 2a confidence levels for the Gold+CMB+BAO dataset using fiom = 0.28 ± 0.03. The upper left hand panel shows 
the confidence levels in Ai — A2, with the black dot representing ACDM. The upper right hand panel shows the logarithmic 
2cr variation of the DE density in terms of redshift. The dashed line represents ACDM. The lower left and right hand panels 
represent the variation of the equation of state and deceleration parameter respectively. The dashed lines in both panels 
represent ACDM. The thick solid line in the lower right hand panel shows the acceleration epoch, i.e. the redshift at which the 
universe started accelerating. 
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FIG. 7: 2a confidence levels for the SNLS+CMB+BAO dataset using fi 0m = 0.28 ± 0.03. The upper left hand panel shows 
the confidence levels in Ai — A2, with the black dot representing ACDM. The upper right hand panel shows the logarithmic 
2cr variation of the DE density in terms of redshift. The dashed line represents ACDM. The lower left and right hand panels 
represent the variation of the equation of state and deceleration parameter respectively. The dashed lines in both panels 
represent ACDM. The thick solid line in the lower right hand panel shows the acceleration epoch, i.e. the redshift at which the 
universe started accelerating. Results are shown upto redshift z = 1.01 



